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Abstract 

The non-photonic neutrino exchange mechanism of the lepton flavor violating 
muon-electron conversion in nuclei is revisited. First we determine the nu¬ 
cleon coupling constants for the neutrino exchange mechanism in a relativistic 
quark model taking into account quark confinement and chiral symmetry re¬ 
quirements. This includes a new, previously overlooked tree-level contribution 
from neutrino exchange between two quarks in the same nucleon. Then for 
the case of an additional sterile neutrino we reconsider the coherent mode of 
this process. The presence of a mixed sterile-active neutrino state zy, heavier 
than the quark confinement scale A c ~ 1 GeV may significantly improve the 
prospects for observation of this process in future experiments as compared to 
the conventional scenario with only light neutrinos. Turning the arguments 
around we derive new experimental constraints on zy* — mixing from the 
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non-observation of muon-electron conversion. 

The recent results from the Super-Kamiokande JTJ and SNO |2j] experiments on atmo¬ 
spheric and solar neutrinos give convincing evidence on neutrino oscillation and hence on 
neutrino masses and lepton flavor violation (LFV). As is known, experimental searches for 
rare processes offer complimentary information on the LFV. This can shed additional light 
on the physics underlying this phenomenon, discriminating various models beyond the stan¬ 
dard model (SM) (for a review see [|3|]). The muon-electron [(/i _ ,e _ )] conversion in nuclei 

M, 

l^b + (A Z) (A, Z) + e , (1) 


is one of the most prominent lepton flavor changing reaction. Experiments searching for 
this process in the coherent channel with a monoenergetic final state electron have reached 
an unprecedented level of sensitivity. Presently, the most stringent upper limits on the 
branching ratio R /ie related to the conversion has been set by the SINDRUM II 

collaboration 
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( 2 ) 


where r /ie and are the rates of the (/r~,e“) conversion and ordinary muon capture, 
respectively. Future experiments will significantly improve these limits. There are proposals 
of the SINDRUM II collaboration to reduce the current limits on the ratio R^ e for 48 Ti 
and 197 Au down to 10 -14 and 6 x 1CV 13 ||, respectively. A new Muon Electron Conversion 
(MECO) experiment on 27 Al is planned at BNL m with an expected sensitivity on the 
branching ratio of about 2 x 10 -17 . Another future project PRIME Jll[ for the (/i~,e _ ) 
conversion on 48 Ti is going to reach a sensitivity of 1CU 18 . The realization of these projects 
would allow to set new stringent constraints on the LFV interactions relevant for the (/i~, e~) 
conversion. This process can be triggered by the LFV interactions associated with the 
exchange of neutrinos and/or new heavy particles (neutralinos, charginos, leptoquarks etc.) 
predicted in models beyond the SM f^-0. In general all the (/r~, e ~) conversion mechanisms 
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can be separated into photonic and non-photonic ones. A mechanism is photonic if it involves 
a virtual photon line connecting the effective leptonic LFV current with the electromagnetic 
nuclear current, otherwise a mechanism is non-photonic. These classes of mechanisms differ 
significantly on their particle and nuclear physics sides and are usually studied independently. 

In this letter we concentrate on the non-photonic neutrino exchange mechanism. We are 
studying a model with three left-handed, weak doublet neutrinos v' Li = (v' Le , v' L/i , v' Lr ) and a 
certain number n of the SM singlet, right-handed sterile neutrinos v Ri = (v R1 , ...v' Rn ). Due 
to mixing they form n + 3 neutrino mass eigenstates N t with masses m; related to the weak 
eigenstates v' a = (v' Le , v' LfX , v' L t , v R1 ,..., v Rn ) by an unitary mixing matrix U as N % — U* a v' a . 

Among Ni there must be at least three observable light neutrinos dominated by the 
active v' e T components while the other states may be of arbitrary mass. In particular, 
they may include additional light neutrinos (zy), one of which might be relevant for the 
phenomenology of neutrino oscillations , as well as intermediate and heavy mass neutrinos 
( v h ). The presence or absence of these neutrino states is an issue for experimental searches. 
For simplicity we consider the (zy, zyj scenario with an arbitrary number of light neutrinos 
zy with masses rn m on the eV scale and one neutrino state vh with mass rrih larger than the 
typical hadronic scale A c ~ 1 GeV. A similar scenario, but with zy, having a mass below 1 
GeV was previously studied in connection with K ± and r semileptonic decays 

The analysis of the nuclear (p~, e~) conversion starts with the elementary nucleon process 
/i~ + N —> A^ + e~. In models with non-trivial neutrino mixing this process can be realized 
at the quark level according to the diagrams of Fig. |l]. The diagrams of Fig. la,b are the 
well-known one-quark box diagrams || while that of Fig. lc is the new tree-level two-quark 
diagram. These lowest order diagrams represent the complete set of the neutrino exchange 
diagrams on quark-level relevant for the above nucleon process. 

In the process considered the typical momentum transfer Q 2 to the nucleon is small 
comparable to the scale set by the muon mass with Q 2 ~ m jp Therefore, the quarks in 
the diagrams considered cannot be treated as free particles as would be the case in the 
asymptotic region Q 2 A^. An appropriate treatment should deal with quarks as states 
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which are confined in the nucleon. Due to the lack of a rigorous theory for confinement in 
QCD one has to engage phenomenological models. In this work we are using the perturbative 
chiral quark model (PCQM) |[hi| . [T4|| treating quarks as extended objects, the constituent 
quarks, which are confined in the nucleon. In this model each quark vertex acquires a 
form factor with the characteristic momentum scale A c ~ 1 GeV related to the confinement 
length l c ~ A" 1 . In the diagrams of Fig. 1 these form factors set the scale for the loop 
momentum q u of the virtual neutrino. This is in contrast to the previous analysis || of 
the diagrams Fig. la,b associated with different (neutrino mass dependent) terms of the 
(/r - ,e - ) conversion amplitude, where the q v scale is set by the W-boson mass due to the 
presence of the corresponding propagators in these diagrams. 

Knowing the characteristic scale q 0 ~ A c of the neutrino momentum q v in the diagrams 
of Fig. 1, we consider the general structure of the (/i _ ,e _ ) conversion amplitude A fie . In 
the (zy, Uh) neutrino scenario introduced above one can write 


A^ e ~ 


E 


77 . 77 *. 


q 2 — m 2 + ie 

^ light 


■ G(q /%)d q 


( 3 ) 


/ rn f u r 

— E U I* U * (! + E - ••• ) ' G (9 2 /?o) A-/ G {q 2 /ql)d 4 q ( m h > qo)- 

q i q rn h j 

Here G{q 2 /q^) is a characteristic function suppressing the contribution for q 2 q% in the 
loop momentum integration. Then it follows that 


(S, C, i (7;,^)/4G( 9 2 / 9o 2 )*, for m h ^q„ 

Al Me ~ < V ° J (4) 

{ -U „ h u* eh $ G(q 2 /ql)^ for m h > q 0 

as a consequence of the unitarity of the mixing matrix with = 0. 

Previously, only the case nih <C qo of Eq. (^) was considered in the literature |||5|, 
Because of the smallness of the ratio m 2 / q q the neutrino exchange mechanism leads to rates 
for the (n ~, e~) conversion which are out of reach for ongoing and near future experiments. 
The situation changes if there exists a heavy neutrino state Vh with mass nih qo and 
with a non-vanishing admixture of active flavors v^ e . In this case the suppression factor 
associated with the small neutrino masses is replaced by the product of mixing matrix 
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elements U^,hU* h as indicated in Eq. (f|). Since the existing experimental constraints on 
LJfj,h,U* h are not stringent |16| one may expect much larger rates for (p _ ,e _ ) conversion in 
the (yi, Vh) scenario than for the case without an intermediate mass i/h state. 

Following the standard approach [0] we consider the effective nucleon Lagrangian written 


as 


£ e J e f '(x) = G 2 f ml U^ h U* h e(x)7 a (l - * 

[POb a (/v - f V Als)v{x) +n( x h a (fv - fAlb)n{x)} + h.c.. 


( 5 ) 


Here is the mass of the muon. The partial contributions of the diagrams Fig. la,b,c are 
contained in the coupling coefficients as fy A = fy A a + fy ' A + fv,A C - I 11 th e present work 
we restrict ourselves to the dominant mode of (p - ,e - ) conversion, where the axial-vector 
current contribution JT5| is neglected and, thus, only the vector form factors fy ,n are relevant 
for the subsequent analysis. 

We evaluated the form factors fy n within the perturbative chiral quark model (PCQM), 

for the study of low-energy 


a relativistic quark model suggested in [13] and extended in 


properties of baryons. The model operates with relativistic quark wave functions and takes 
into account quark confinement as well as chiral symmetry requirements. The PCQM was 
successfully applied to cr-term physics and to the electromagnetic properties of the nucleon 
|T4|1 . In the present analysis we included the contributions from both the one-body (Fig. 


UK b) and the two-body (Fig. [Tjc) diagrams neglecting the external three-momenta of the 
leptons. For the one-body diagrams of Fig. 0a and Fig. [lb we restrict the expansion of the 
quark propagator to the ground state eigenmode: 


iGA*, y) -» iG 0 (x, y) = u 0 (x)u 0 (y)e ,a,( " K) 0(x o -y 0 ), 


( 6 ) 


where S 0 and u 0 (x) are the quark ground state energy and wave function; that is we restrict 
the intermediate baryon states to N and A configurations. In Ref. |T|] we showed that this 
approximation for the quark propagator works quite well in the phenomenology of low-energy 
nucleon physics. 
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With above approximations the partial contributions of the diagrams of Fig. la,b,c to 


the coupling constant of the vector current are: 


nN;la _ - 

JV ~ o 


d 3 k 

(2^)3 


d 3 x / d 3 y 


f N;lb _ 
JV ~ 


\k\ — mn — ie 

3 / _ ^ \ (*) 
x < N\J2{ $o(+)$o(+) + $(+)$(+)) \N >, 
»= 1 ' 3 

d 3 k 


( 7 ) 


(2tt) 


d 3 x / d 3 y e ?;fc( ^7) 


|fc| + 

(*) 


( 8 ) 


fAf;lc _ 

Jv — 


x < W|^(<F 0 (f)<h 0 (+)+ $(+)$(+)) |AT >, 

Z— 1 

d 3 k 


f AJL f d z x f d 3 y ( , ---- ) 

J (27t) 3 J J V \k\ — m,. — ie \k\ + m,.J 


\k\ — rriu — ie \k\ + 


(9) 


x < W|^($o(x) w $o(x) (i) + $(f) (?;) $(f) ( ^j|iV> . 

<=# V ' 

Here f and j are the quark indices, <f>o(+) = uo(x)Xo u oY) and &(x) = uo(x)juo(x) are 
the time and spatial components of the quark vector current. The single components are 
projected onto the three-quark state building up the nucleon state \N >. As in Refs. |T4| 
we use the variational Gaussian ansatz [T7|] for the quark ground state wave function given 
by: 


( 


Uq(x) = N exp 


x‘ 


2 R 2 


\ 


in HfL 

Y p R 


Xs Xf Xc : 


( 10 ) 


where Xs, Xf: Xc refer to the spin, flavor and color spinors. The constant N = [7r 3 / 2 f? 3 (l + 
3p 2 /2 )]~ 1//2 is fixed by the normalization condition. 

Our Gaussian ansatz contains two model parameters: the dimensional parameter R 
and the dimensionless parameter p. The parameter p can be related to the axial coupling 
constant (Ja calculated in zeroth-order (or 3q-core) approximation: 


9 a = - 1 


2 P 2 


1 + 1 P 2 


51 + 2y 
3 3 ’ 


( 11 ) 


where 7 is a relativistic reduction factor 


1 - 9 1 

7 = , L = —9a ~ -■ 


1 + 1 P 1 
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( 12 ) 
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The parameter R can be physically understood as the mean radius of the three-quark core 
and is related to the charge radius < r% > EO °f the proton in the leading-order (or zerotli- 
order) approximation as [[14] 

3 R 2 1 + f p 2 


< r E >LO — 


2 i + |p 2 "" V" 2 


= R 2 { 2 


(13) 


In our calculations we use the tree-level value g\— 1-25 as obtained in Chiral Perturbation 
Theory |U| and the averaged value of R = 0.6 fin |14| corresponding to < r E > EO = 0.6 fm 2 . 

A straightforward analytical evaluation of the expressions in Eqs. (|7|)-(|9]) results in the 
following values for the partial isospin dependent vector coupling constants: 


f$ la = 2.37 + i 0.41, 
ff b = 0.64, 


fy c = 0.44 + i 0.14, /y ;ic = 0.44 + i 0.14, 

For the total nucleon vector coupling constants fy' n , entering in Eq. (||), we finally obtain 
the values 


fy’ la = 1.19 + i 0.21, 
fv ]U = 1 - 27 , 


n;lc 


(14) 


fP = 3.45 + i 0.55, fy = 2.90 + 7 0.35, 


(15) 


which are just the sum of the partial contributions. 

Starting from the effective Lagrangian of Eq. (0) the branching ratio of the coherent 
(p~,e~) conversion is derived as 

2 


R^e = K h u:„\ 2 - m „ ( G F m ‘) 


2\4 Pe-^e jy( ry ^ \ || 


F(Z,p e ) 


/ jt ^ m 2 \ 71 l'' 

n in fL 1 


(16) 


Here, E e (E e = — £&, Eb is the muon binding energy) and p e (p e = \p e \) are energy and 


momentum of the outgoing electron. F{Z,p e ) is the relativistic Coulomb factor [19] and the 
nuclear structure factor is defined as 

1 




(ftM p + fZM r 


(17) 


In our analysis we used values for the nuclear matrix elements M. P) n as derived in Ref. 
@. Using the calculated vector coupling constants as an input, in Table |T| we indicate the 
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numerical values for the relevant quantities entering in Eq. (|16]). With these values we get 
the following results for the coherent (/r _ ,e _ ) conversion on the nuclear targets 4H Ti, 197 Au 
and 27 Al as: 


R 


theor 

fie 


= 3.19 x Kn n ( 2 M(), 6.34 x l(r 11 ( 48 T'i), 4.73 x l(T 10 ( 187 .4!i), 


|£V4*J 2 


(18) 


which is valid for a heavy neutrino state zy, with mass nth 3> A c ~ 1 GeV. From these results 
we derive upper limits on the product of mixing matrix elements which correspond 

to the sensitivity of present and near future experiments discussed in the introduction. These 
limits which are set by the experimental upper bounds are listed in Table ||. The present 
constraint \U tl hU* h \ < 0.1, provided by the SINDRUM 11 48 Ti experiment, is rather weak. An 
improvement on this bound is expected from the ongoing SINDRUM II 197 Au experiment. 
Significant sensitivity, down to 10 -4 — 10 -3 , will be hopefully achieved by the future MECO 
(target 27 Al) and PRIME (target 48 Ti ) experiments. 

The mixing of massive neutrinos, like Vh with active v e ,n,T flavors was previously looked 
for in various experiments except for (fi~, e ~) conversion. An extensive list for the constraints 
on the U eh ( and | U^l mixing matrix elements for various masses nth is given in Ref. |TB 


For nth < 19-6 GeV all the values of U lt h and | U e h have been excluded by the MARK II 
collaboration. The ALEPH collaboration ruled out all values of these matrix elements for 
25.0 < nth A 42.7 GeV. However, above 42.7 GeV there exist only narrow domains of nth 
where \U a h\, with a = e,/i, are constrained. The typical constraints are \U Q h\ 2 < 10~ 10 — 
10 -5 . The (/i _ , e ~) conversion constraints in Table [TJ cover the partially constrained region 
nth > 42.7 GeV and extend into the region nth > 1600 GeV currently unconstrained. On 
the other hand, the poorly explored region nth > 42.7 GeV offers a loop-hole for observation 
of the (/r~, e ~) conversion in future experiments as its rate is unconstrained in the most part 
of this mass region. 

In summary, we have studied the non-photonic neutrino exchange mechanism of coherent 
conversion in nuclei in the presence of sterile neutrinos. We found a new tree level 
contribution to the (/i~,e _ ) conversion (Fig. lc) which is as important as the previously 






known box-type contributions (Fig. la,b). The nucleon form factors, parameterizing the 
effective nucleon Lagrangian, have been analyzed within the perturbative chiral quark model 
II . In this model the momentum scale of the virtual neutrino is set by the quark confine¬ 
ment with A c ~ 1 GeV in the three diagrams. This significantly differs from the previous 
analysis [3] of the diagrams Fig. la,b where this scale is of the order of ~ M\y. The lowering 
the neutrino momentum scale has a notable effect on the analysis of the observability of the 
(p~, e_ ) conversion. We have shown that in the neutrino scenario with at least one heavy 
neutrino state Uh with mass A c ~ 1 GeV the rate of this lepton flavor violating process 

could be large enough to be observed in planed experiments. This observation is in contrast 
to the conventional belief that (/i _ ,e _ ) conversion will not be detected experimentally even 
in the distant future if the process is dominated by the neutrino mechanism. It turn we 
derived new upper bounds on the product of neutrino mixing matrix elements from the 
non-observation of this process in running and planed experiments. 
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TABLES 

TABLE I. The nuclear structure factor (dimensionless) and other quantities of Eq. (|l6|). 
The experimental values of the total muon capture rate r /ti , are from Ref. [20]. 


Nucleus 

Pe 

(MeV) 

iWtxio 6 *- 1 ) 

F(Z,p e ) 

|A4 M e| 

27 A1 

0.531 

-0.470 

0.71 

1.37 

0.754 

48 Ti 

0.529 

-1.264 

2.60 

1.62 

1.87 

197 Au 

0.485 

-9.938 

13.07 

3.91 

7.37 


TABLE II. Upper bounds on the product of mixing matrix elements | U /t hU* h | of a heavy 

neutrino with v e ,, flavors in the mass region rrih S> 1 GeV as derived from the sensitivity of 

present and near future conversion experiments. 


Present limits 


Expected limits 

Nucleus 

R™r- Ref. \U, h U* eh \ 

R™p- 

Ref. 

K h u* eh \ 

27 Al 


< 5 x nr 17 

|10] 

< 1.2 x 10~ 3 

48 Ti 

< 6.1 x 10~ 13 § < 9.8 x 10" 2 

< 1 x 10” 18 

111] 

< 1.3 x 10~ 4 

197 Au 

< 2.0 x 10- 11 § < 0.21 

< 6 x 10~ 13 

i 

< 3.5 x 10 -2 
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FIGURES 



FIG. 1. Feynman diagrams for (/i , e ) conversion in nuclei associated with neutrino exchange 
on the quark level. Both one- (a,b) and two-body (c) mechanisms are considered. 
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